Pathogenic effects of Leu200Pro and Arg387His VRK1 protein variants on phosphorylation targets and H4K16 acetylation in distal hereditary motor neuropathy.

Rare recessive variants in the human VRK1 gene are associated with several motor neuron diseases (MND), such as amyotrophic lateral sclerosis, spinal muscular atrophy, or distal hereditary motor neuropathies (dHMN). A case with dHMN carrying two novel VRK1 gene variants, expressing Leu200Pro (L200P) and Arg387His (R387H) variant proteins, identified that these protein variants are functionally different. The Leu200Pro variant shares with several variants in the catalytic domain the loss of the kinase activity on different substrates, such as histones, p53, or coilin. However, the distal Arg387His variant and the distal Trp375* (W375X) chinese variant, both located at the end of the low complexity C-terminal region and proximal to the termination codon, retain their catalytic activity on some substrates, and mechanistically their functional impairment is different. The L200P variant, as well as most VRK1 pathogenic variants, impairs the phosphorylation of BAF and histone H4K16 acetylation, which are required for DNA attachment to the nuclear envelope and chromatin accessibility to DNA repair mechanisms, respectively. The R387H variant impairs phosphorylation of H2AX, an early step in different types of DNA damage responses. The functional variability of VRK1 protein variants and their different combinations are a likely contributor to the clinical phenotypic heterogeneity of motor neuron and neurological diseases associated with rare VRK1 pathogenic variants. KEY MESSAGES: VRK1 variants implicated in motor neuron diseases are functionally different. The L200P variant is kinase inactive, and the R387H variant is partially active. VRK1 variants alter H4K16 acetylation and loss of coilin and BAF phosphorylation. VRK1 variants alter Cajal bodies and DNA damage responses. VRK1 variant combination determines the neurological phenotype heterogeneity.

STAG2: Computational Analysis of Missense Variants Involved in Disease.

The human STAG2 protein is an essential component of the cohesin complex involved in cellular processes of gene expression, DNA repair, and genomic integrity. Somatic mutations in the STAG2 sequence have been associated with various types of cancer, while congenital variants have been linked to developmental disorders such as Mullegama-Klein-Martinez syndrome, X-linked holoprosencephaly-13, and Cornelia de Lange syndrome. In the cohesin complex, the direct interaction of STAG2 with DNA and with NIPBL, RAD21, and CTCF proteins has been described. The function of STAG2 within the complex is still unknown, but it is related to its DNA binding capacity and is modulated by its binding to the other three proteins. Every missense variant described for STAG2 is located in regions involved in one of these interactions. In the present work, we model the structure of 12 missense variants described for STAG2, as well as two other variants of NIPBl and two of RAD21 located at STAG2 interaction zone, and then analyze their behavior through molecular dynamic simulations, comparing them with the same simulation of the wild-type protein. This will allow the effects of variants to be rationalized at the atomic level and provide clues as to how STAG2 functions in the cohesin complex.

Computational Modeling and Design of New Inhibitors of Carbapenemases: A Discussion from the EPIC Alliance Network.

The EPIC consortium brings together experts from a wide range of fields that include clinical, molecular and basic microbiology, infectious diseases, computational biology and chemistry, drug discovery and design, bioinformatics, biochemistry, biophysics, pharmacology, toxicology, veterinary sciences, environmental sciences, and epidemiology. The main question to be answered by the EPIC alliance is the following: "What is the best approach for data mining on carbapenemase inhibitors and how to translate this data into experiments?" From this forum, we propose that the scientific community think up new strategies to be followed for the discovery of new carbapenemase inhibitors, so that this process is efficient and capable of providing results in the shortest possible time and within acceptable time and economic costs.

Diversity of mechanisms to control bacterial GTP homeostasis by the mutually exclusive binding of adenine and guanine nucleotides to IMP dehydrogenase.

IMP dehydrogenase(IMPDH) is an essential enzyme that catalyzes the rate-limiting step in the guanine nucleotide pathway. In eukaryotic cells, GTP binding to the regulatory domain allosterically controls the activity of IMPDH by a mechanism that is fine-tuned by post-translational modifications and enzyme polymerization. Nonetheless, the mechanisms of regulation of IMPDH in bacterial cells remain unclear. Using biochemical, structural, and evolutionary analyses, we demonstrate that, in most bacterial phyla, (p)ppGpp compete with ATP to allosterically modulate IMPDH activity by binding to a, previously unrecognized, conserved high affinity pocket within the regulatory domain. This pocket was lost during the evolution of Proteobacteria, making their IMPDHs insensitive to these alarmones. Instead, most proteobacterial IMPDHs evolved to be directly modulated by the balance between ATP and GTP that compete for the same allosteric binding site. Altogether, we demonstrate that the activity of bacterial IMPDHs is allosterically modulated by a universally conserved nucleotide-controlled conformational switch that has divergently evolved to adapt to the specific particularities of each organism. These results reconcile the reported data on the crosstalk between (p)ppGpp signaling and the guanine nucleotide biosynthetic pathway and reinforce the essential role of IMPDH allosteric regulation on bacterial GTP homeostasis.

The clinical and molecular spectrum of QRICH1 associated neurodevelopmental disorder.

De novo variants in QRICH1 (Glutamine-rich protein 1) has recently been reported in 11 individuals with intellectual disability (ID). The function of QRICH1 is largely unknown but it is likely to play a key role in the unfolded response of endoplasmic reticulum stress through transcriptional control of proteostasis. In this study, we present 27 additional individuals and delineate the clinical and molecular spectrum of the individuals (n = 38) with QRICH1 variants. The main clinical features were mild to moderate developmental delay/ID (71%), nonspecific facial dysmorphism (92%) and hypotonia (39%). Additional findings included poor weight gain (29%), short stature (29%), autism spectrum disorder (29%), seizures (24%) and scoliosis (18%). Minor structural brain abnormalities were reported in 52% of the individuals with brain imaging. Truncating or splice variants were found in 28 individuals and 10 had missense variants. Four variants were inherited from mildly affected parents. This study confirms that heterozygous QRICH1 variants cause a neurodevelopmental disorder including short stature and expands the phenotypic spectrum to include poor weight gain, scoliosis, hypotonia, minor structural brain anomalies, and seizures. Inherited variants from mildly affected parents are reported for the first time, suggesting variable expressivity.

Dysfunctional Homozygous VRK1-D263G Variant Impairs the Assembly of Cajal Bodies and DNA Damage Response in Hereditary Spastic Paraplegia.

BACKGROUND AND OBJECTIVES: To conduct a genetic and molecular functional study of a family with members affected of hereditary spastic paraplegia (HSP) of unknown origin and carrying a novel pathogenic vaccinia-related kinase 1 (VRK1) variant. METHODS: Whole-exome sequencing was performed in 2 patients, and their parents diagnosed with HSP. The novel VRK1 variant was detected by whole-exome sequencing, molecularly modeled and biochemically characterized in kinase assays. Functionally, we studied the role of this VRK1 variant in DNA damage response and its effect on the assembly of Cajal bodies (CBs). RESULTS: We have identified a very rare homozygous variant VRK1-D263G with a neurologic phenotype associated with HSP and moderate intellectual disability. The molecular modeling of this VRK1 variant protein predicted an alteration in the folding of a loop that interferes with the access to the kinase catalytic site. The VRK1-D263G variant is kinase inactive and does not phosphorylate histones H2AX and H3, transcription factors activating transcription factor 2 and p53, coilin needed for assembly of CBs, and p53 binding protein 1, a DNA repair protein. Functionally, this VRK1 variant protein impairs CB formation and the DNA damage response. DISCUSSION: This report expands the neurologic spectrum of neuromotor syndromes associated with a new and rare VRK1 variant, representing a novel pathogenic participant in complicated HSP and demonstrates that CBs and the DNA damage response are impaired in these patients.

Clinical relevance of postzygotic mosaicism in Cornelia de Lange syndrome and purifying selection of NIPBL variants in blood.

Postzygotic mosaicism (PZM) in NIPBL is a strong source of causality for Cornelia de Lange syndrome (CdLS) that can have major clinical implications. Here, we further delineate the role of somatic mosaicism in CdLS by describing a series of 11 unreported patients with mosaic disease-causing variants in NIPBL and performing a retrospective cohort study from a Spanish CdLS diagnostic center. By reviewing the literature and combining our findings with previously published data, we demonstrate a negative selection against somatic deleterious NIPBL variants in blood. Furthermore, the analysis of all reported cases indicates an unusual high prevalence of mosaicism in CdLS, occurring in 13.1% of patients with a positive molecular diagnosis. It is worth noting that most of the affected individuals with mosaicism have a clinical phenotype at least as severe as those with constitutive pathogenic variants. However, the type of genetic change does not vary between germline and somatic events and, even in the presence of mosaicism, missense substitutions are located preferentially within the HEAT repeat domain of NIPBL. In conclusion, the high prevalence of mosaicism in CdLS as well as the disparity in tissue distribution provide a novel orientation for the clinical management and genetic counselling of families.

DLG4-related synaptopathy: a new rare brain disorder.

PURPOSE: Postsynaptic density protein-95 (PSD-95), encoded by DLG4, regulates excitatory synaptic function in the brain. Here we present the clinical and genetic features of 53 patients (42 previously unpublished) with DLG4 variants. METHODS: The clinical and genetic information were collected through GeneMatcher collaboration. All the individuals were investigated by local clinicians and the gene variants were identified by clinical exome/genome sequencing. RESULTS: The clinical picture was predominated by early onset global developmental delay, intellectual disability, autism spectrum disorder, and attention deficit-hyperactivity disorder, all of which point to a brain disorder. Marfanoid habitus, which was previously suggested to be a characteristic feature of DLG4-related phenotypes, was found in only nine individuals and despite some overlapping features, a distinct facial dysmorphism could not be established. Of the 45 different DLG4 variants, 39 were predicted to lead to loss of protein function and the majority occurred de novo (four with unknown origin). The six missense variants identified were suggested to lead to structural or functional changes by protein modeling studies. CONCLUSION: The present study shows that clinical manifestations associated with DLG4 overlap with those found in other neurodevelopmental disorders of synaptic dysfunction; thus, we designate this group of disorders as DLG4-related synaptopathy.

Pathogenic convergence of CNVs in genes functionally associated to a severe neuromotor developmental delay syndrome.

BACKGROUND: Complex developmental encephalopathy syndromes might be the consequence of unknown genetic alterations that are likely to contribute to the full neurological phenotype as a consequence of pathogenic gene combinations. METHODS: To identify the additional genetic contribution to the neurological phenotype, we studied as a test case a boy, with a KCNQ2 exon-7 partial duplication, by single-nucleotide polymorphism (SNP) microarray to detect copy-number variations (CNVs). RESULTS: The proband presented a cerebral palsy like syndrome with a severe motor and developmental encephalopathy. The SNP array analysis detected in the proband several de novo CNVs, nine partial gene losses (LRRC55, PCDH9, NALCN, RYR3, ELAVL2, CDH13, ATP1A2, SLC17A5, ANO3), and two partial gene duplications (PCDH19, EFNA5). The biological functions of these genes are associated with ion channels such as calcium, chloride, sodium, and potassium with several membrane proteins implicated in neural cell-cell interactions, synaptic transmission, and axon guidance. Pathogenically, these functions can be associated to cerebral palsy, seizures, dystonia, epileptic crisis, and motor neuron dysfunction, all present in the patient. CONCLUSIONS: Severe motor and developmental encephalopathy syndromes of unknown origin can be the result of a phenotypic convergence by combination of several genetic alterations in genes whose physiological function contributes to the neurological pathogenic mechanism.

Novel Dominant KCNQ2 Exon 7 Partial In-Frame Duplication in a Complex Epileptic and Neurodevelopmental Delay Syndrome.

Complex neurodevelopmental syndromes frequently have an unknown etiology, in which genetic factors play a pathogenic role. This study utilizes whole-exome sequencing (WES) to examine four members of a family with a son presenting, since birth, with epileptic-like crises, combined with cerebral palsy, severe neuromotor and developmental delay, dystonic tetraparexia, axonal motor affectation, and hyper-excitability of unknown origin. The WES study detected within the patient a de novo heterozygous in-frame duplication of thirty-six nucleotides within exon 7 of the human KCNQ2 gene. This insertion duplicates the first twelve amino acids of the calmodulin binding site I. Molecular dynamics simulations of this KCNQ2 peptide duplication, modelled on the 3D structure of the KCNQ2 protein, suggest that the duplication may lead to the dysregulation of calcium inhibition of this protein function.

VRK1 (Y213H) homozygous mutant impairs Cajal bodies in a hereditary case of distal motor neuropathy.

BACKGROUND: Distal motor neuropathies with a genetic origin have a heterogeneous clinical presentation with overlapping features affecting distal nerves and including spinal muscular atrophies and amyotrophic lateral sclerosis. This indicates that their genetic background is heterogeneous. PATIENT AND METHODS: In this work, we have identified and characterized the genetic and molecular base of a patient with a distal sensorimotor neuropathy of unknown origin. For this study, we performed whole-exome sequencing, molecular modelling, cloning and expression of mutant gene, and biochemical and cell biology analysis of the mutant protein. RESULTS: A novel homozygous recessive mutation in the human VRK1 gene, coding for a chromatin kinase, causing a substitution (c.637T > C; p.Tyr213His) in exon 8, was detected in a patient presenting since childhood a progressive distal sensorimotor neuropathy and spinal muscular atrophy syndrome, with normal intellectual development. Molecular modelling predicted this mutant VRK1 has altered the kinase activation loop by disrupting its interaction with the C-terminal regulatory region. The p.Y213H mutant protein has a reduced kinase activity with different substrates, including histones H3 and H2AX, proteins involved in DNA damage responses, such as p53 and 53BP1, and coilin, the scaffold for Cajal bodies. The mutant VRK1(Y213H) protein is unable to rescue the formation of Cajal bodies assembled on coilin, in the absence of wild-type VRK1. CONCLUSION: The VRK1(Y213H) mutant protein alters the activation loop, impairs the kinase activity of VRK1 causing a functional insufficiency that impairs the formation of Cajal bodies assembled on coilin, a protein that regulates SMN1 and Cajal body formation.

Delineation of phenotypes and genotypes related to cohesin structural protein RAD21.

RAD21 encodes a key component of the cohesin complex, and variants in RAD21 have been associated with Cornelia de Lange Syndrome (CdLS). Limited information on phenotypes attributable to RAD21 variants and genotype-phenotype relationships is currently published. We gathered a series of 49 individuals from 33 families with RAD21 alterations [24 different intragenic sequence variants (2 recurrent), 7 unique microdeletions], including 24 hitherto unpublished cases. We evaluated consequences of 12 intragenic variants by protein modelling and molecular dynamic studies. Full clinical information was available for 29 individuals. Their phenotype is an attenuated CdLS phenotype compared to that caused by variants in NIPBL or SMC1A for facial morphology, limb anomalies, and especially for cognition and behavior. In the 20 individuals with limited clinical information, additional phenotypes include Mungan syndrome (in patients with biallelic variants) and holoprosencephaly, with or without CdLS characteristics. We describe several additional cases with phenotypes including sclerocornea, in which involvement of the RAD21 variant is uncertain. Variants were frequently familial, and genotype-phenotype analyses demonstrated striking interfamilial and intrafamilial variability. Careful phenotyping is essential in interpreting consequences of RAD21 variants, and protein modeling and dynamics can be helpful in determining pathogenicity. The current study should be helpful when counseling families with a RAD21 variation.

Evaluating Face2Gene as a Tool to Identify Cornelia de Lange Syndrome by Facial Phenotypes.

Characteristic or classic phenotype of Cornelia de Lange syndrome (CdLS) is associated with a recognisable facial pattern. However, the heterogeneity in causal genes and the presence of overlapping syndromes have made it increasingly difficult to diagnose only by clinical features. DeepGestalt technology, and its app Face2Gene, is having a growing impact on the diagnosis and management of genetic diseases by analysing the features of affected individuals. Here, we performed a phenotypic study on a cohort of 49 individuals harbouring causative variants in known CdLS genes in order to evaluate Face2Gene utility and sensitivity in the clinical diagnosis of CdLS. Based on the profile images of patients, a diagnosis of CdLS was within the top five predicted syndromes for 97.9% of our cases and even listed as first prediction for 83.7%. The age of patients did not seem to affect the prediction accuracy, whereas our results indicate a correlation between the clinical score and affected genes. Furthermore, each gene presents a different pattern recognition that may be used to develop new neural networks with the goal of separating different genetic subtypes in CdLS. Overall, we conclude that computer-assisted image analysis based on deep learning could support the clinical diagnosis of CdLS.

MEPSAnd: minimum energy path surface analysis over n-dimensional surfaces.

SUMMARY: n-dimensional energy surfaces are becoming computationally accessible, yet interpreting their information is not straightforward. We present minimum energy path surface analysis over n-dimensional surfaces (MEPSAnd), an open source GUI-based program that natively calculates minimum energy paths across energy surfaces of any number of dimensions. Among other features, MEPSAnd can compute the path through lowest barriers and automatically provide a set of alternative paths. MEPSAnd offers distinct plotting solutions as well as direct python scripting. AVAILABILITY AND IMPLEMENTATION: MEPSAnd is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSAnd/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A novel RAD21 p.(Gln592del) variant expands the clinical description of Cornelia de Lange syndrome type 4 - Review of the literature.

Cornelia de Lange syndrome (CdLS) is a heterogeneous developmental disorder where 70% of clinically diagnosed patients harbor a variant in one of five CdLS associated cohesin proteins. Around 500 variants have been identified to cause CdLS, however only eight different alterations have been identified in the RAD21 gene, encoding the RAD21 cohesin complex component protein that constitute the link between SMC1A and SMC3 within the cohesin ring. We report a 15-month-old boy presenting with developmental delay, distinct CdLS-like facial features, gastrointestinal reflux in early infancy, testis retention, prominent digit pads and diaphragmatic hernia. Exome sequencing revealed a novel RAD21 variant, c.1774_1776del, p.(Gln592del), suggestive of CdLS type 4. Segregation analysis of the two healthy parents confirmed the variant as de novo and bioinformatic analysis predicted the variant as disease-causing. Assessment by in silico structural model predicted that the p.Gln592del variant results in a discontinued contact between RAD21-Lys591 and the SMC1A residues Glu1191 and Glu1192, causing changes in the RAD21-SMC1A interface. In conclusion, we report a patient that expands the clinical description of CdLS type 4 and presents with a novel RAD21 p.(Glu592del) variant that causes a disturbed RAD21-SMC1A interface according to in silco structural modeling.

Human Mitochondrial HMG-CoA Synthase Deficiency: Role of Enzyme Dimerization Surface and Characterization of Three New Patients.

Mitochondrial 3-hydroxy-3-methylglutaryl-CoA synthase deficiency (mitochondrial HMG-CoA synthase deficiency or mHS deficiency, OMIM #605911) is an inborn error of metabolism that affects ketone body synthesis. Acute episodes include vomiting, lethargy, hepatomegaly, hypoglycemia and dicarboxylic aciduria. The diagnosis is difficult due to the relatively unspecific clinical and biochemical presentation, and fewer than 30 patients have been described. This work describes three new patients with mHS deficiency and two missense mutations c.334C>T (p.R112W) and c.430G>T (p.V144L) previously not reported. We developed a new method to express and measure the activity of the enzyme and in this work the study is extended to ten new missense variants including those of our patients. Enzymatic assays showed that three of the mutant proteins retained some but seven completely lacked activity. The identification of a patient homozygous for a mutation that retains 70% of enzyme activity opens the door to a new interpretation of the disease by demonstrating that a modest impairment of enzyme function can actually produce symptoms. This is also the first study employing molecular dynamics modelling of the enzyme mutations. We show that the correct maintenance of the dimerization surface is crucial for retaining the structure of the active center and therefore the activity of the enzyme.

A nucleotide-controlled conformational switch modulates the activity of eukaryotic IMP dehydrogenases.

Inosine-5'-monophosphate dehydrogenase (IMPDH) is an essential enzyme for nucleotide metabolism and cell proliferation. Despite IMPDH is the target of drugs with antiviral, immunosuppressive and antitumor activities, its physiological mechanisms of regulation remain largely unknown. Using the enzyme from the industrial fungus Ashbya gossypii, we demonstrate that the binding of adenine and guanine nucleotides to the canonical nucleotide binding sites of the regulatory Bateman domain induces different enzyme conformations with significantly distinct catalytic activities. Thereby, the comparison of their high-resolution structures defines the mechanistic and structural details of a nucleotide-controlled conformational switch that allosterically modulates the catalytic activity of eukaryotic IMPDHs. Remarkably, retinopathy-associated mutations lie within the mechanical hinges of the conformational change, highlighting its physiological relevance. Our results expand the mechanistic repertoire of Bateman domains and pave the road to new approaches targeting IMPDHs.

Two-step ATP-driven opening of cohesin head.

The cohesin ring is a protein complex composed of four core subunits: Smc1A, Smc3, Rad21 and Stag1/2. It is involved in chromosome segregation, DNA repair, chromatin organization and transcription regulation. Opening of the ring occurs at the "head" structure, formed of the ATPase domains of Smc1A and Smc3 and Rad21. We investigate the mechanisms of the cohesin ring opening using techniques of free molecular dynamics (MD), steered MD and quantum mechanics/molecular mechanics MD (QM/MM MD). The study allows the thorough analysis of the opening events at the atomic scale: i) ATP hydrolysis at the Smc1A site, evaluating the role of the carboxy-terminal domain of Rad21 in the process; ii) the activation of the Smc3 site potentially mediated by the movement of specific amino acids; and iii) opening of the head domains after the two ATP hydrolysis events. Our study suggests that the cohesin ring opening is triggered by a sequential activation of the ATP sites in which ATP hydrolysis at the Smc1A site induces ATPase activity at the Smc3 site. Our analysis also provides an explanation for the effect of pathogenic variants related to cohesinopathies and cancer.

A Practical Quantum Mechanics Molecular Mechanics Method for the Dynamical Study of Reactions in Biomolecules.

Quantum mechanics/molecular mechanics (QM/MM) methods are excellent tools for the modeling of biomolecular reactions. Recently, we have implemented a new QM/MM method (Fireball/Amber), which combines an efficient density functional theory method (Fireball) and a well-recognized molecular dynamics package (Amber), offering an excellent balance between accuracy and sampling capabilities. Here, we present a detailed explanation of the Fireball method and Fireball/Amber implementation. We also discuss how this tool can be used to analyze reactions in biomolecules using steered molecular dynamics simulations. The potential of this approach is shown by the analysis of a reaction catalyzed by the enzyme triose-phosphate isomerase (TIM). The conformational space and energetic landscape for this reaction are analyzed without a priori assumptions about the protonation states of the different residues during the reaction. The results offer a detailed description of the reaction and reveal some new features of the catalytic mechanism. In particular, we find a new reaction mechanism that is characterized by the intramolecular proton transfer from O1 to O2 and the simultaneous proton transfer from Glu 165 to C2.